Finite temperature scaling theory for the collapse of Bose-Einstein condensate 
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We show how to apply the scahng theory in an inhomogeneous system Uke harmonically trapped 
Bose condensate at finite temperature. We calculate the temperature dependence of the critical 
number of particles by a scaling theory within the Hartree-Fock approximation and find that there 
is a dramatic increase in the critical number of particles as the condensation point is approached. 
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I. INTRODUCTION 

Experimentalpl i, i, i, i, i, 0] and theoretical [1 
B Eyj HH, studies of Bose-Einstein condensa- 

tion(BEC) reveal the fact that Bose-Einstein conden- 
sates of trapped ultracold atomic gases are influenced by 
atomic interaction 

B 1 [M [13, [H [3 which in general 
is characterized by s-wave scattering length(as) and can 
be tuned [H, [H, [o, [ll] by an external magnetic field. 
Stability and collapse of the Bose gas with negative scat- 
tering length has been observed in the clouds of ultracold 
■^Lifo,!!!! and s'^Rbfli.llOl. If the interaction is attrac- 
tive, the gas tends to increase the density of the central 
region of the trap. This tendency is opposed by the zero- 
point energy and thermal energy of the atoms. If the 
number of atoms is greater than a critical numbcr(A*'c), 
the central density increases strongly and the zero-point 
and thermal energies are no longer able to avoid the col- 
lapse of the gas. 

Let us first consider a many particle system of 3-D 
isotropic harmonically trapped ideal Bose gas of N alkali 
atoms. The system is in equilibrium with its surround- 
ings at temperature (T). Let the mass of each particle be 
m. The ground state of a single particle is macroscop- 
ically populated below a certain temperature which is 
called the Bose-Einstein condensation temperature and 
is given by Q Tc = ^[-^V^^, where oj is the angu- 
lar frequency of the harmonic oscillators and k is the 
Boltzmann constant. Now we consider the gas to be in- 
teracting weakly and attractively by the two body in- 
teracting potential V^mt(r) = gS^{r), where g = — "^^^ ° 
is the coupling constant and a = —as is the absolute 
value of the s-wave scattering length. For < T < Tc 

which 



the length scale of the system becomes I — 
is the length scale of the ground state of the harmonic 
oscillators. The typical two body interaction energy for 
N number of particles is[3l ~ ^^5/2^^, where V is the 
volume of the system. For the attractive interaction, the 
gas tends to increase the density of the central region of 



the trap. For T = 0, this tendency is resisted by the 
zero-point energy of the atoms. If the number of atoms 
is greater than a critical number, the central density in- 
creases so much that the zero-point energy is no longer 
able to avoid the collapse of the gas. In the critical situa- 
tion, the typical oscillator energy must be comparable to 
the typical interaction energy and as a consequence we 
have ^ - 1. 

However, for < T < Tc, although length scale of the 
system is / yet the thermal fluctuation dominates over the 
quantum fluctuation. In this situation, the collapse of the 
attractive Bose gas is opposed by the thermal energy (Q 
3?iw[-|^]'*C(4) hwN^/^) of the atoms and consequently, 
we get the estimation of the critical number as ~ 

[111/2 > ^ 

However, close to the condensation temperature, the 
length scale of the system becomes Q Lt ~ ^\f^- For 
this reason the estimation of the critical number would 
be ^ - [{f > 1. 

In this paper we will explicitly calculate the critical 
numbers by a scaling theory within the Hartree-Fock(H- 
F) approximation and will explicitly show its tempera- 
ture dependence. Before doing that, let us see how the 
critical number for T = 0, was calculated by the scaling 
theory of Baym and Pethick • 



II. ZERO TEMPERATURE SCALING THEORY 
FOR THE COLLAPSE OF EEC 



For T ^ 0, the system of N weakly interacting Bose 
particles is well described by the Gross-Pitaevskii(G-P) 
equation 



ih 



5*o(r, 



dr 



(- 



2m 



V{v) 
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I vl/„(r,r) nvI/o(r,r), 
(1) 

where ^(r) = ^moPr"^ is the harmonic potential and 
^o(r,T) is the order parameter for Bose condensate and 
is a function of space(r) and time(T). A stationary so- 
lution of the G-P equation have the form ^'□(r, r) — 
4'o(r)e~T^ where /i is the chemical potential. The sta- 
tionary solution of this form can be obtained by mini- 



2 



mizing the energy functional Q 

= /(|^ I VvI/o(r) p +Vir) I vI;o(r) p +| | vI/o(r) \'^)d'r 

(2) 

for a fixed number of particles (A^ = / I P dr). The 
stationary solution without current gives ^'o(r) — \/n(r), 
where n{r) is the number density of particles at the po- 
sition r. For stationary and for no current condition, we 
have the energy functional from the Eqn.(2) as 



E{n) = 




For the noninteracting gas(g = 0), the solution of G-P 
equation must be the ground state of harmonic oscilla- 
tors. Hence, for g = 0, the solution of the G-P equation 

is \J 13^3/2 ■ For <? < 0, the solution of G-P equation 

was obtained by a scaling theory (variational approach) 
based on the Gaussian function [lol. [T^ . From this scal- 
ing theory we can write *o(i") = ^fn = y^^^j^^e^^^ , 
where ^' is a scaling parameter which fixes the width of 
the condensate. Since the attractive interaction causes 
the reduction of the width of the trap, v must be less 
than one. With this solution, the energy functional of 
the Eqn.(3) can be written in units of Nhuj as 

X^{v) = {E(y)lNh^\ = + v^) - 4^^]. (4) 

According to FIG. 1, this energy functional has a 
minimum and a maximum below a critical number of 
particles(Afc)- In equilibrium, the system must be at the 
minimum of the energy functional. About this minimum, 
V has a stable value. Above the critical number of parti- 
cles, E(y) has no stable minimum. In this situation, the 
energy of the system arbitrarily reduces to take = 
and the gas collapses to vanish its width. Nc is calcu- 
lated by requiring that the first and second derivative 
vanish at the critical point. From this scaling theory we 
can have d [TO, [il ^ =0.671. According to the numer- 
ical result [2]| = 0.575 and according to the experi- 
mental result[2^ — 0.459. However, the sarne prob- 
lem was also attacked by many authors [12, [l^, [3 [2^ . 
Kagan, Surkov and Shlyapnikov discussed the dynamic 
properties of a trapped Bose-condensed gas under vari- 
ations of the confining field and found analytical scal- 
ing solutions for the evolving condensate 22]. Pitaevskii 
addressed the condition of collapse from the dynami- 
cal point of viewf23'|. Yukalov and Yukalova calculated 
the optimal trap shape allowing for the condensation 
of the maximal number of atoms with negative scatter- 
ing lengths[23|. Savage, Robins and Hope solved the 
Gross-Pitaevskii equation numerically for the collapse, 
induced by a switch from positive to negative scattering 
lengths [25|. Although the G-P equation is extensively 
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FIG. 1: Plot of energy per particle {Xo{u))in units hcu with 
the scaling parameter v. For the thick line, = 0.671. For 
the dotted line, ^ = 0.5. For the dashed line, ^ = 0.35. 
All the lines follow from the Eqn.(4). 



used for the Bose gas with attractive interaction, yet it is 
to be noted that, it was rigorously derived for repulsive 
interaction[26j. 



III. SCALING ASSUMPTION FOR THE 
COLLAPSE OF BEC AT FINITE TEMPERATURE 

Although the problem of an attractive Bose gas at fi- 
nite temperature was discussed by many authors [27l [28l. 
[29| . yet the temperature dependence of critical number 
of particles with a scaling theory has not been calculated. 
The authors of the Refs.(27l. [28j considered the collapse 
of the condensate only. The authors of the Ref. i^Sj ex- 
plored the collapse for T = and for T ^ Tc. But we 
are considering the collapse for <T <Tc and are con- 
sidering the fact that both the condensed cloud and the 
thermal cloud of the Bose gas would collapse due to the 
attractive type of interaction. In the fifth section we will 
calculate Nc by the following scaling theory within the 
Hartree-Fock(H-F) approximation which was discussed 
extensively for Bose gas in the Ref.0. 

For the interacting Bose gas, average occupation num- 
bers fii as well as single particle wave functions are 
determined [g| by minimizing the grand potential Q = 
{E -TS - iiN), where 5* = fc^^(l + n,)Zn(l -t- n,) - 
{ni)ln{ni) is the entropy and i runs from to oo. The 
entropy(S') and the total number of particles(A^ — rii) 
both are the functions of {rii}. Minimizing the grand 
potential with respect to the occupation number (rii) we 
get the average occupation number as Ui = ^(g._^f//;T — j- 



3 



Since the entropy and the total number of particles de- 
pend on the occupation numbers {n^}, the single par- 
ticle wave functions {(pi} are simply obtained by mini- 
mizing the energy functional with proper normalization 
constraint / d^r | P= 1 d]. Within the Hart ree-Fock 
approximation, we have the expression of energy func- 
tional asQ 



#0 



+V{r)no{r) + V(r)nT(r) + ^nl(r) 
+2gno{r)nTir) + gn^r)] 



(5) 



Minimizing the above energy functional of the Eqn.(5) 
we have the following equations [1] . 



2m 



[-^V2 + Vir) + 2gn(r)]0, 



(7) 



From the Eqns.(6) and (7) we can determine the wave 
functions and the energy eigenvalues {e^} for a given 
no(r) and riTir). In our scaling theory we introduce 
a scaling parameters z/, which would fix the width of 
the condensed and thermal clouds as well as would take 
proper choices of no(r) and nT(r)- Eventually the choices 
of the wave functions are denoted by the scaling pa- 
rameter v. Since the minimum of the grand potential 
corresponds to the equilibrium of the system, the mini- 
mum of the energy functional for a certain choice of {4>i} 
would correspond to the equilibrium of the system. So, 
for equilibrium condition of the system, the energy func- 
tional would be minimized with respect to the scaling 
parametcr(i^). In the FIG. 2, we will see that above a 
critical number of particles(A^c), the energy functional 
has no stable minimum. If the energy functional as well 
as the grand potential has no stable minimum, the grand 
potential as well as the energy functional would arbitrar- 
ily decrease until the width of the system becomes zero. 
Under this condition the system is said to be collapsed. 

In absence of the interaction, the condensed cloud is 
described by the wave function (boM = 

and by the density no(r) = 13^3/2 ^' ■ In ab- 
sence of interaction, the number density of the excited 
particles isQ ^^(r) = "^^^ ''^ — 



r + - 



(27r?i)3 



.3/33/2 



), where At = yj^: and 53(0;) = x + 
. is Bose-Einstein function of a real 



1^9 lie' 
xV23/2 + a;" 

variable x. The length scale of the excited particles in 

Since 



thermal equilibrium is Lt ~ y = I x 
the length scale(LT) of the thermal cloud is proportional 



to the length scale(Z = y/h/moj) of the condensed cloud, 
the two length scale would reduce by the same factor due 
to the attractive interaction. 

In presence of the attractive interaction the wave func- 
tion of the condensed would become 



Mr) 



1 



and the density of the thermal cloud would become 

1 



1/3 A3 



(8) 



(9) 



The same scaled form of nrir) is also obtained from 
the variation of B-E statistics such that n(p, r) — 
— . With this variational form of statis- 



,2^2 



+ V{r) + g{no{r) + 2nT(r))]</)o = ficj^o (6) Nt = J 



tics we have the total number of excited particles as 
(e-f^) = (|^)3C(3). ForT^T,, 



there is no particles in the ground state and we have 
N = (7^)3C(3). For T < T^, the total number of parti- 
cles in the ground state is rig = iV — Nt = N{1 — 3^)^. 



IV. EVALUATION OF THE HARTREE-FOCK 
ENERGY FUNCTIONAL 

Let us now evaluate the seven terms of energy func- 
tional in the Eqn.(5). Since we have scaled the length 
and there is a square of the derivative of the wave func- 
tions with respect to the length, all the kinetic energy 
terms of Eqn.(5) would be proportional to the inverse 
squared of the scaling parameter. Due to the presence 
of the term quadratic in length in the expression of har- 
monic potential energy, all the potential energy terms of 
Eqn.(5) would be proportional to the squared of the scal- 
ing parameter. Similarly, as the interaction energy goes 
like the inverse of the volume, all the interacting energy 
terms of Eqn.(5) would be proportional to the inverse 
cubed of the scaling parameter. 

The first term of the energy functional is the kinetic 
energy of the condensed cloud and is obtained as 



no 47r 

2^I/3;37r3/2 (,.2/2)2 
4 J/2 Tn 



2 -r^/u^l^ 2j 

r e ' ' r dr 



(10) 



The second term of the energy functional is the kinetic 
energy of the thermal cloud and is approximately ob- 
tained as 



00 ^00 ^2 



1 



4Vrfp, 2 



JO '^rn^(^+^^)/kT _^{2^hy- 



Anrr dr 



(11) 
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The third term of the energy functional is the potential 
energy of the condensed cloud and is obtained as 



collecting all the terms of the energy functional from the 
Eqn.(lO) to (16), we have 



1 



HQ 



(12) 



The fourth term of the energy functional is the potential 
energy of the thermal cloud and is obtained as 



1 



J 2 2 

[^1 (e~ 2u'^kT )]-mLV^r^4:7Tr^dr 



r (c(3) 



2^ ^ ^4/3/^^4 



(13) 



The fifth term of the energy functional is the interaction 
energy of the particles solely in the condensed cloud and 
is obtained as 



m (1/3/3^3/2)2 



Airr dr 



(14) 



The sixth term of the energy functional is the interac- 
tion energy of the particles in the condensed and thermal 
clouds and is obtained as 



no 



1 



m 1^3/3^^3/2 j/3^3 



Z/2 



23/^VC(372) a,,3/2(J^)3/2[^ 



(f)']^- (15) 



The seventh term of the energy functional is the inter- 
action energy of the particles solely in the thermal cloud 
and is obtained as 



E 



TO ly^X^ ;n 



j3/2j3/2 



dr 



2 S'N^/^ a 7^9/2 



(16) 



where S' = EZ=i (^j)V-^l,+j)V2 « 0.6534. 

The energy functional(ii') of the Eqn.(5) is now a func- 
tion of the variational parameter v. Above the critical 
number of particles, E{y) has the lowest value(— oo) at 
V = Q. So, above the critical number of particles, the 
width of the system would be zero to signalize the col- 
lapse. Let us define the energy per particle in units of 
hw as Xt{v) — Let us also write ^ = t. Now 



"8C(3/2)jVV^a ^3/,^^ ^3 




am 



(17) 



where 

and 

C2(i) = 



3 N^/^Ci4^)t* 



(18) 



1 Na 



+ S' 




8C(3/2) TVV^a 3/2. 



(19) 



From the Eqns.(18) and (29), and putting f = 0.0066 
[20 l| and t — 0.8, in the Eqn.(17) we can write 



Xo.sM = (0.366-l-0.52037Vi/3)(^ -l-i^2) 



0.0069iVi/2 + 0.000627iV 



(20) 



The Eqn.(20) represents the FIG. 2. In the FIG. 2 we 
see that below a critical number of particles, Xtiv) has 
a stable minimum and an unstable maximum. At the 
critical number of particles, the minimum and maximum 
coincide. Above the critical number of particles, there is 
no stable minimum and the energy of the system arbi- 
trarily decreases to — cxi and its width(i/^) becomes zero 
to signalize the collapse of the gas. 



V. CRITICAL CONDITIONS FOR THE 
COLLAPSE 

For a minimum of Xtiv) as well as for the equilibrium 
condition of the system, we must have = 0. At the 
critical point(t/c), the minimum and maximum coincide. 
So, at the critical point, we must have |i/^= 0. For 

the critical point(i'c), the Eqn.(17) gives following equa- 
tions 

f)X 9 

^ I..- -ci(i)-3 + 2ci(t)z. + c,{t)-^ = (21) 
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FIG. 2: Plot of energy per particle (Xo.8(/^))in units huj with 
the scaling parameter u. For the thick line, = 27.78. For 



I 

the dashed line, = 20. For the dotted line, 
the lines follow from the Eqn.(20). 
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FIG. 3: Plot of the critical number A'^c in units of - with 

a 

temperature(r) in units of (Tc). This plot follows form the 
Eqn.(25). 



and 



\u=ci {t) A + 2ci (t) - C2 (t) = 0. (22) 



From the Eqns.(21) and (22) we have i/c = (1/5)^/** = 
0.66874 and C2{t) = ^i'cCi{t), which represents the equa- 
tion for the critical number(A^c) and can be recast as 




312 



8C(3/2)iVy'a 



3/2 r 



7r[C(3)]3 I ~ 15 =U^^ ^ 



3 7Vy^C(4)i^ 
'2 [C(3)]4/3 



(23) 



From the Eqn.(23) we get the expression of for 
t = as 



Nca 
I 



0.671 



(24) 



Now we see our scaling theory for T > to be consisted 
with that given by Baym and Pethick for T = [l^ . 

For < t < 1, we can consider the second and third 
terms of the left hand side and the first term of the right 
hand side of the Eqn.(23) as the perturbing terms. We 
expand its left hand side with respect to its first term 
and similarly, we expand its right hand side with respect 
to its second term. Now we get the critical number from 



the Eq.(23) up to three leading orders in Z/a as 



Nca 
I 



^^Milll!^! [1^13/2 rli 1/2 

[C(3)]2 ^ 5 ^ ^a^ [l-t3]3 
4[l-<3] 



3[ /8C(^,3/2(i_,3)^^, 



mm 



^t9/2 



[27r]V«[C(4)]3/4 4Z. 3^ g ^3 
[C(3)] ^ 5 ^ [l-i3]V2 



1.210[i-]i/2^ 



1.006 



a' (1-^3)3 [l„i3] 



- [10.666i3/2 



r [i-t3]V2 



(25) 



The Eqn.(25) represents the behavior of the critical num- 
ber with temperature. ForfSo*! f = 0.0066, the change of 
critical number(A'c) with t is shown in FIG. 3. From this 
figure and according to the Eqn.(25), we get ^ = 27.89 
for t = 0.8. The Eqn.(25) represent an approximate re- 
lation between and t. However, from the FIG. 2, we 
see that ^ = 27.78 for t = 0.8. 

At the condensation temperature i.e. as t —f 1, only 
the third term of the left hand side of the Eqn.(23) con- 
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tributes to give 

/ ^flJ ^2^ S"6 ^ 5 ^ 

= 2.253[-]^ (26) 
a 

From the above equation we see that near the condensa- 

1 /6 

tion temperature j(= 1.145) is of the order of unity 
and is also expected from the quahtative point of view. 

VI. CONCLUSIONS 

Initially we explained the physics of the collapse of 
the attractive atomic Bose gas. Then we qualitatively 
estimated Nca/l for the various range of temperatures 
(0 < T < Tc). Finally we calculate Nca/l by a scal- 
ing theory within the Hartree-Fock approximation. Our 
calculation supports the qualitative estimation of Nca/l. 

In the FIG. 3, we see that the temperature dependence 
of Nc is not significant for < T < 0.5Tc. If experiment 
of collapse is performed within this range of tempera- 
ture, the zero temperature theory would well satisfy the 
experimental result. 

Here we consider an attractive type of delta potential 
which is not an actual potential but a model potential 
■ The problem with this model potential is that, 
no bound state is formed for a finite number of particles 
in a three dimensional delta potential. So, for a true delta 
potential, the collapse of a finite number of particles is 
impossible. However, once we do the scaling assumption, 
the collapse is not impossible for a delta potential and we 
can calculate the critical numbers. 

It is a coincidence of our scaling theory that the grand 
potential and the energy are minimized simultaneously 



for a certain choice of the scaling parameter. However, 
in general, the grand potential and the energy are not 
minimized simultaneously. If we want to calculate the 
critical number by some other theory, the minimization 
of the energy may not be sufRcient for the calculation 
of the critical number. In that case, some modification 
would come from the entropy and chemical potential. 

Although the Hartree-Fock approximations are exten- 
sively usedjl, even for the critical region(Tc) of Bose 
systems, yet the Hartree-Fock-Bogoliubov(H-F-B) ap- 
proximation is a more reliable one ^2^]. Any mean- field 
approximation(H-F, H-F-B, etc.) is also known to be less 
suitable for the critical region. 

We considered the interaction as Virit(r) = g5'^{Y). For 
repulsive interaction, the density of the gas in the central 
region decreases. This would cause a slight decrease in 
Tc as shown theoretically in the Ref.(l3| and experimen- 
tally in the Refs. 0, [s^l . On the other hand, for attractive 
interaction, the density in the central region of the trap 
would increase and we expect slight increase in Tc- How- 
ever, our scaling theory does not predict any Tc shift. 
This is a limitation of our scaling theory. We should look 
for a better theory which would simultaneously predict 
the critical numbers and the Tc shift. 
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